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Abstract. PHENIX preliminary results on the J/ 1 ^ production in Cu + Cu and 
Au + Au collisions at y/sjvw = 200 GeV are presented. They are compared to results 
from lower energy experiments NA50 and NA60 at CERN SPS and to expectations 
from various theoretical models. 
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O . 1. Introduction 

"ft 

Based on the assumption that cc pairs are created in primordial hard nucleon-nucleon 
collisions only (due to the large mass of the charm quark), J/\l/'s suppression by color 
screening in a deconfmed medium [1] is suggested to be used as a signature of a Quark- 
Gluon Plasma (QGP) formation in the relativistic heavy ion collisions. However, at 
RHIC energies, up to N cS ~ 10-20 cc pairs can be created in the most central Au + Au 
collisions. Since the probability to recombine uncorrelated c and c goes as N cE , in- 
medium formation of J/\l/'s is foreseen at RHIC, leading to an enhancement of the J/\I/ 
yield [2, 3]. Alternate mechanisms, not involving QGP, such as the interaction with 
hadronic co- movers [4] (which are secondaries produced in the collision), also successfully 
accounted for the abnormal J/^> suppression in most central Pb + Pb collisions, in excess 
of the normal nuclear absorption, as seen by the lower energy experiment NA50 [5, 6] at 
CERN SPS. The measurements at higher energy could help disentangling this variety of 
models. In addition to these effects which can be expected in case of the QGP formation, 
other so-called cold nuclear effects [7] are known to affect the J/\l/ yield, in particular: 
(i) the shadowing (initial-state effect) which is a modification of the parton distribution 
function of a free nucleon by the nuclear environment, and (ii) the nuclear absorption 
(final-state effect) where the formed cc pre-resonance is broken by interactions with 
primary target /projectile nucleons. d + Au measurements is used to evaluate the cold 
nuclear effects at RHIC. 
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Figure 1. Left panel: J /if? nuclear modification factor measured in d + Au as 
a function of rapidity and compared to theoretical expectations for cold nuclear 
effects [8]. Kopeliovich [9] and FGS parametrization use a stronger gluon shadowing 
than the EKS98 parametrization [10, 11]. Right panel: J/^f? nuclear modification 
factor measured in d + Au, Cu + Cu and Au + Au vs the number of participating 
nucleons [12] compared to Vogt [7] predictions in Au + Au collisions for cold nuclear 
effects. Vertical bars stand for statistical errors, brackets for point-to-point and boxes 
for global systematic errors. 



2. Results and discussion 



During the RHIC data taking periods Run-3 to Run-5 (2003-2005) PHENIX measured 
the J/* yield in p + p, d + Au, Cu + Cu and Au + Au at ^/s]^ = 200 GeV. J/* 
measurements are performed in two decay channels, J/^> — > e + e~ at mid rapidity 
(\y\ < 0.35), and J /\& — > yU + yU~ at forward rapidity (1.2 < \y\ < 2.2) intervals. Here, we 
will focus on Cu + Cu and Au + Au preliminary results (resp. Run-5 and Run-4). Run-3 
p+p results [8] define the baseline of the J/ty production and hence are used to compute 
the nuclear modification factor defined as Rab = M^dN^/dy wnere N co ii is the number 
of binary collisions. Rab is expected to be unity if there is no effect. The Run-3 d + Au 
results [8] are used to characterize the cold nuclear effects at RHIC. The left panel of 
figure 1 shows the nuclear modification factor RdAu as a function of rapidity. At positive 
rapidity (<i-going direction) are probed Au partons carrying a small momemtum fraction 
x, which is the region where the gluon shadowing is expected. RdAu is indeed significantly 
lower than unity at positive rapidity. Theoretical expectations [9, 10, 11] for nuclear 
absorption and shadowing effects are shown for comparison. Data favors a modest 
shadowing in agreement with the EKS98 parametrization and a modest absorption 
cross-section a a b s = 1 mb, but a stronger absorption (3 mb) is still compatible within 
our large errors. A higher luminosity d + Au data set is desirable to better constrain 
the cold nuclear effects at RHIC. The preliminary results for the centrality dependence 
of the nuclear modification factor in Cu + Cu and Au + Au [12] are shown on the 
right panel of figure 1, together with the published minimum bias d + Au value. The 
suppression pattern are compared to the predictions [7] which use EKS98 shadowing 
parametrization for central (y = 0) and forward rapidities (y = 2) and assume 1 mb and 
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3 mb nuclear absorption. Even relative to the 3 mb nuclear absorption case, the J/^> 
suppression in most central Au + Au collisions goes significantly beyond the cold nuclear 
effects at RHIC. A factor « 3 of suppression is seen relative to binary scaled p + p. This 
is surprisingly close to the suppression factor seen at lower -s/snn = 17.3 GeV by the 
NA50 experiment [5, 6] for the most central Pb + Pb collisions relative to the NA51 p+p 
result [13]. Moreover, suppression at RHIC energy follows the same general trend [14] as 
NA50 (and consequently preliminary NA60 [15]) suppression, in spite of very different 
conditions, namely s/^nn, cold nuclear effects (at SPS, a abs = 4.18 ± 0.35 mb), and 
the estimated maximum energy density (~ 3 GeV/fm 3 at SPS, 5 GeV/fm 3 at RHIC, 
if assuming the same formation time of 1 fm/c). On the left panel of figure 2, the 
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Figure 2. J '/^f nuclear modification factor in d + Au, Cu + Cu and Au + Au vs 
the number of participating nucleons [12] compared to various models of final-state 
interaction in the medium. Left panel: solid line stands for the co- mover model [4], 
dashed lines for a direct suppression in a hot medium [2]. Right panel: solid lines is 
the sum of the direct suppression (dashed lines) and the regeneration [2] (dotted lines). 

centrality dependence of the nuclear modification factor is compared to two different 

models of final-state interactions that successfully described the NA50 results: the first 

one suppresses J/Ws using co-movers interactions [1], and the second involves direct 

dissociation in a hot medium [2]. Both models seem to overestimate the J/^> suppression 

at RHIC energy. On the right panel of figure 2 the charm quark recombination 

mechanism [2] is added to the direct dissociation and the sum is in better agreement 

with the observed suppression. To explore further the recombination scenario, more 

observables are compared to the related predictions from [3] on figure 3, namely the 

centrality dependence of the mean square transverse momentum < p\ > (left and middle 

panel) and the evolution with centrality of the J/ \1/ yield as a function of rapidity (right 

panel). In this model, recombined J/ty are predicted to have a lower < p\ > since charm 

quark production is supposed to be significant at low p?- Observed values of < p\ > 

have large errors and lay between the two hypotheses, with and without recombination. 

J/ty's from recombination should also have a rapidity spectra peaked at y = where 

the charm density is expected to be larger. This implies a narrowing rapidity spectra 

with an increasing centrality, which seems not being observed here. 
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Figure 3. J mean square transverse momentum as a function the number of binary 
collisions in p + p, d + Au, Cu + Cu (left panel) and Au + Au [12] (middle panel). 
It is compared to Thews predictions [3] with recombination (solid line) and without 
recombination (dashed line). Open makers are for the mid rapidity measurements, and 
solid markers for the forward rapidity ones. Right panel: J/vI 7 yield as a function of 
rapidity and increasing centrality from p + p to Cu + Cu to most central Au + Au [12]. 
Vertical bars stand for statistical errors and boxes for point-to-point systematic errors. 

3. Summary 

PHENIX preliminary results on J/ty production measured at mid and forward rapidities 
in Cu + Cu and Au + Au at y/s^N = 200 GeV show a suppression pattern well 
beyond cold nuclear effects, but still similar to the suppression seen at CERN SPS. 
Models that successfully describe SPS results overestimate the suppression at RHIC. 
Adding recombination seems to restore the agreement for the centrality dependence, but 
conclusions from testing recombination against over observables (mean square transverse 
momentum, shape of the rapidity spectra) remain unclear. 
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